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Abstract 

We have studied the finite A^-unit Langevin model subjected to multiplicative 
noises, by using the augmented moment method (AMM), as a continuation of our 
previous paper [H. Hasegawa, J. Phys. Soc. Jpn. 75 (2006) 033001]. Effects of cou- 
plings on stationary and dynamical properties of the model have been investigated. 
The difference and similarity between the results of diffusive and sigmoid couplings 
are studied in details. Time dependences of average and fluctuations in local and 
global variables calculated by the AMM are in good agreement with those of di- 
rect simulations (DSs) . We also discuss stationary distributions of local and global 
variables with the use of the Fokker-Planck equation (FPE) method and DSs. It is 
demonstrated that stationary distributions show much variety when multiplicative 
noise and external inputs are taken into account. 
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1 INTRODUCTION 



The Langevin equation has been widely employed as a useful model for a wide range of 
stochastic phenomena. Much study has been made on the Langevin model for a single unit 
as well as coupled systems (for a recent review, see Ref. [1]). The Langevin equation has 
been commonly solved by using the Fokker-Planck equation (FPE) method For A^-unit 
Langevin equations, the FPE method leads to (A^ + l)-dimensional partial equations to be 
solved with proper boundary conditions, which is usually very difficult. Direct simulation 
(DS) requires the computational time which grows as A^^ with increasing N. As a useful 
semi-analytical method for stochastic equations, Rodriguez and Tuckwell j3] proposed 
the moment method in which the first and second moments of variables are taken into 
account. In this approach, original A^- dimensional Langevin equations are transformed 
to {N/2){N + 3)-dimensional deterministic equations. For example, this figure becomes 
65 and 5150 for = 10 and A^ = 100, respectively. Based on a macroscopic point of 
view, Hasegawa 4J has proposed the augmented moment method (AMM), in which the 
dynamics of coupled Langevin equations is described by a fairly small number (three) of 
quantities: averages and fiuctuations of local and global variables. The AMM has been 
successfully applied to a study on the dynamics of coupled stochastic systems described 
by Langevin, FitzHugh-Nagumo and Hodgkin-Huxley models subjected to additive noises 
with global, local or small-world couplings (with and without transmission delays) 
The AMM was originally developed by expanding variables around their mean values in 
the stochastic model in order to obtain the second-order moments both for local and global 
variables In a recent paper PU]; we have reformulated the AMM with the use of the 
FPE, in order to apply the AMM to coupled Langevin model subjected to multiplicative 
noises, in which the difficulty of the Ito versus Stratonovich representations is inherent. 

In recent years, much attention has been paid to multiplicative noises in addition 
to additive noises (for a review of study on multiplicative noises, see Ref. related 
references therein). The stationary distribution of the Langevin model subjected to mul- 
tiplicative noises has been considerably investigated in various contexts Inter- 
esting phenomena caused by the two noises have been intensively studied. It has been 
realized that the property of multiplicative noises is different from that of additive noises 
in some respects. (1) Multiplicative noises induce the phase transition, creating an ordered 
state, while additive noises are against the ordering [in]-|2ni- (2) Although the proba- 
bility distribution in stochastic systems subjected to additive Gaussian noise follows the 
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Gaussian, multiplicative Gaussian noises generally yield non-Gaussian distribution jl2j- 
(3) The scaling relation of the effective strength for additive noise given by 
P{N) = I3{1)/\/N is not applicable to that for multiplicative noise: a{N) ^ a{l)/\fN, 
where a{N) and (3{N) denote effective strengths of multiplicative and additive noises, 
respectively, in the A^-unit system A naive approximation of the scaling relation for 
multiplicative noise: a{N) = a{V)/\fN as adopted in Ref. [20] yields the result which 
does not agree with that of DSs. 

The purpose of the present paper is to discuss effects of couplings on stationary and 
dynamical properties of the A^-unit Langevin model with multiplicative noises, which 
has been not investigated in jTUI. The paper is organized as follows. In Section 2, the 
AMM is employed for a discussion on the finite- A^ Langevin model which is subjected to 
additive and multiplicative noises and which is coupled by diffusive and sigmoid couplings. 
Numerical results are presented in Section 3. Section 4 is devoted to discussion and 
conclusion, where the stationary distribution of local and global variables are studied 
with the FPE and DS. 

2 AUGMENTED MOMENT METHOD 

2.1 A Generalized Langevin model 

We have adopted the finite A^-unit Langevin model given by 

^ = + aG{x,)r^,{t) + PUt) + l['\t) + &\t), (1) 

at 



with 



and 



I^if) = 7 E k'c(t) - xm + y E H{x,(t)), {^ = l-N) (2) 



H{x) = -4=. (3) 

Here F{x) and G{x) denote arbitrary functions of x: J and K express the diffusive and 
sigmoid couplings, respectively, whose effects will be separately discussed in Sections 2.2 
and 2.3: Z (= N — 1) stands for the coordination number: I^^\t) is an external input: 
a and (3 denote the strengths of multiplicative and additive noises, respectively, and rjiit) 
and ^i{t) express zero-mean Gaussian white noises with correlations given by 

{v^{t)vAt')) = s^At't'), (4) 
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imW)) = S,,6{t-t'), (5) 

{mit)Ut')) = 0. (6) 

Although various types of sigmoid functions such as tanh(x) and 1/[1 + exp(x)], etc. have 
been employed in the literature, we here adopt a simple analytical expression given by 
Eq. (3) 

The Fokker-Planck equation for the distribution of p{{xi},t) is given by [21] 

dt ^ dxk 2 

+ lT.f2{[^"G{x,y + [3']p{{x,},t)}, (7) 

where Ik = 1^^ + G'{x) = dG{x)/dx, and = 1 and in the Stratonovich and Ito 
representations, respectively. The averaged, global variable X{t) is given by 

X{t) = j-Y.^^it), (8) 

i 

for which the Fokker-Planck equation P{X, t) is formally given by 

P(X, t)=J---JUidx, p{{x,}, t)6{X-j^Y.^^)- (9) 

i 

We will discuss the property of the coupled Langevin model with the use of the AMM, 
which is the second-moment theory for local and global variables jH^j. The moments 
of local and global variables are defined by 

(xf) = JUidxip{{xi},t)x^i, (10) 

{X'') = j dXP{X,t)X^. (A; = 1,2,-) (11) 

From Eqs. (1), (7), (8), (10) and (11), equations of motions for mean, variance and 
covariance of local variable (xj) and global variable (X) are given by ^01 

^ = (F(a:,)) + (/,) + ^(G'(a:,)G(x,)), (12) 



dt 



-- {xi F{xj)) + {xj F{xi)) + {xilj) + {xjli) 

+ ^[{x:,G\x,)G{x^)) + {x,G\xi)G{xi))] 

+ [a^{G{x,f)+P^]5,,, (13) 

d{X) ^ 1 ^ djxj) 

dt dt ' ^ ' 

d{X^) _ 1 ^^ d{xi Xj) 

dt - dt ' 



where = 1^^^ + /*^^^ Equation (12) is adopted in the mean-field approximation [T7j. In 
Ref. Eqs. (12) and (13) are employed for a discussion on the fiuctuation-induced 

phase transition in infinite-A^ stochastic systems. Equations (14) and (15) play a crucial 
role in discussing finite- iV systems, as will be shown shortly. 

In the AMM jHCn], we take into account the three quantities: /x, 7 and p expressing 
the mean of X, the averaged fluctuations in local variables (xi) and fluctuations in global 
variable (X), respectively, which are deflned by 

p = W = ^E(^^)' (16) 

i 

i 

P = {{X-pf). (18) 

Expanding Xi in Eqs. (12)-(15) around the average value of /x as 

Xi = n + 6xi, (19) 

and retaining up to the order of {6xi6xj), we get equations of motions for /i, 7 and p given 
by 

^ = fo + f2l + K[ho + h2^]+(^^ygogi + 3{gig2 + 9093)1] + (20) 

+ i(P+l){gl + 2gog2)a^l + a^go+P^ (21) 
^ = 2f,p + 2Khp+{<P + l)igl + 2gog2)a'p+^ + ^, (22) 

where fe = {l/£\)d^F{p)/dx\ ge = {l/i\)d^G{p)/dx^, and he = {l/i\)d^H{p)/dx^. Orig- 
inal A^- dimensional stochastic equations given by Eqs. (l)-(3) are transformed to three- 
dimensional deterministic equations given by Eqs. (20)-(22). The stability of Eqs. (20)- 
(22) may be examined by calculating their Jacobian matrix, as will be discussed shortly. 
We note that equations of motions for p and p in Eqs. (20) and (22) do not include the 
term of J for diffusive coupling, while they include the term of K for the sigmoid coupling. 
This is because the average of Y^i{Ii) and in Eqs. (12) and (15) vanish due to the 

nature of the diffusive coupling given in Eq. (2). If we consider the conditional average 
of {xi)i for a given site i, its equation of motion has a term relevant to the coupling J, as 
discussed in Ref. [THj . 
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2.2 Diffusive couplings 

For the linear Langevin model given by F{x) — —Xx and G{x) — x with diffusive cou- 
plings {J ^ 0, K — 0), we get equations of motions for /i, 7 and p in the Stratonovich 
representation (0 = 1) given by 

I = + ^ + (23) 

^ = -2A7 + (^) - 7) + + aV + /3', (24) 

I . _2.,.2.V.^.§. (25) 

The stabihty of the stationary solutions given by Eqs. (23)-(25) may be examined by 
calculating their Jacobian matrix. We get three eigenvalues of A — q;^/2, 2X — 2a'^ + 2JN/Z 
and 2 A — 2q;^, from which the stability condition of the stationary solution is given by 
CK^ < A. The stable stationary solutions for I^^^ — I are given by 

^ = (A-lv2)' ^''^ 

^ («V + /3^)[l + J/Z(A-a^)] 
^ 2(A - q;2 + J7V/Z) ' ^ ' 

- + (asAT^oo) (28) 



2(A -a^ + J) 



2iV(A-a2)' 

p _ ^ / X-a^ + JN/Z \ 
7 ~ iV 1^ X-a^ + J/Z j ' 



(29) 
(30) 



2.3 Sigmoid couplings 

We will make an analysis of the linear Langevin model with F{x) = —Xx and G{x) = x 
for sigmoid couplings {J = 0, K 0) in Eqs. (2) and (3), for which equations of motion 
for II, 7 and p in the Stratonovich representation are given by 

^ = -Xii + ^ + Kiho + h,^) + &\ (32) 

-2A7 + +2q;7 + q;// (33) 



I = -2Ap + 2aV + 2X/iiP+^ + ^, (34) 



where ho = /i/V/i^ + 1, h = + 1)3/2 ^ -(3/i/2)/(/i2 + 1)5/2, r^Yie stationary 

solutions with /'^^^ = I for a small jj, for which H{ii) ~ /x are given by 

^ = (A-aV2-K)' ^^^^ 

{a'^^'' + P^)[l + K/Z{\-a'-K)] 
^ = 2iX -a^ + K/Z) ' 

^ 2(A-a2) ' (-^^-^) (37) 



p 1 / \-a' + K/Z 



7 V^-«^-^ + ^/^ 

1 / X-a^ \ 



(39) 



as iV ^ oo) (40) 



Three eigenvalues of Jacobian matrix relevant to Eqs. (35), (36) and (38) are A — — fC, 
X — — K and A — + K/Z, from which we get the stability condition of < A — i^. 

2.4 Comparison between diffusive and sigmoid couplings 

Comparing Eqs. (23)-(25) for the diffusive coupling with Eqs. (32)-(34) for the sigmoid 
coupling, we note the following difference and similarity in /i, 7 and p. 

(i) /i for the diffusive coupling is independent of the coupling (J) while /i for the sigmoid 
coupling depends on its coupling (K). 

(ii) When the (positive) coupling is introduced, 7 for the diffusive coupling is decreased 
while 7 for the sigmoid coupling is almost independent of it because p ~ 7/A^ for small 
K in the second term of Eq. (33). 

(iii) When the (positive) coupling is introduced, p for the sigmoid coupling is increased 
while p for the diffusive coupling is independent of it. 

(iv) With increasing the (positive) coupling strength, the ratio of p/7 is increased for both 
the couplings. This leads to an increased synchronization ratio of S{t): 

_ ( p{t)h{t) - l/N \ _ ( Np{t)-^{t) \ 

which is one and zero for the completely synchronous and asynchronous states, respec- 
tively 111123. 
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2.5 Nature of the AMM 

Before proceeding to the next section of numerical results, we will discuss the nature of 
the AMM, which is essentially the second-moment approximation for local and global 
variables. One of disadvantages of the AMM is that its applicabihty is limited to the 
weak-noise case because higher-order moments are assumed to be neglected. The second- 
moment given by Eq. (10) is positive definite for the magnitude of the multiplicative 
noise a, as given by 

in the case ol I — J — K — {). A simple calculation leads to the equation of motion of 
the k-th moment for even k given by 

= ^ (x^-^), (A; = 2,4,-) (43) 

from which the stationary value of {x^) is given by 



(A- - 1) 



(45) 



2*^/2 nJ5(A-^Q;2)' 

We get the positive definite {x^) for < 2X/k, which implies that for 2X/k < a"^ < X with 
k > 4, the k-th moment may diverge even if (x^) remains finite. This might throw some 
doubt on the validity of the AMM for the multiplicative noise. Equations (43) expresses 
that the motion of {x'') depends on those of its lower moments (< k — 2), but does not 
on its higher moments (> A; -|- 2). Even if (x^) diverges, for example, it has no effects on 
the motion of (x^) . We hope that our AMM is meaningful and useful for discussions on 
stochastic systems subjected to multiphcative noise, because the results of the AMM are 
in good agreement with those of DS, as will be demonstrated in our numerical calculations. 
The advantage of the AMM is that we can easily discuss the dynamics of A^-unit Langevin 
model by solving the three-dimensional ordinary differential equations. Note that it is 
much more tedious to solve (iV + l)-dimensional partial differential equations in FPE and 
A^-dimensional stochastic equations in DS. 
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3 Numerical results 



3.1 Stationary property 

3.1.1 Size (iV) dependence 

We have performed numerical calculations for linear Langevin models, solving AMM 
equations by the Runge-Kutta method with a time step of 0.01. Direct simulations for 
the A^-unit Langevin model have been performed by using the Heun method with a time 
step of 0.0001. Results shown in the paper are averages of 1000 trials otherwise noticed. 

The A^-dependences of 7 and p in the stationary states for the diffusive couplings 
of J = 0.0, 0.2 and 0.5 are plotted in Fig. 1(a) where sohd curves and marks denote 
the results of the AMM and DS, respectively. We note that for J = 0, p is inversely 
proportional to A^, as realized in Eq. (29). With increasing J, 7 is decreased while p has 
no changes. The results of AMM are in good agreement with those of DS. 

The A^- dependences of 7 and p in the stationary states for sigmoid couplings of A' = 
0.0, 0.2 and 0.5 are plotted in Fig. 1(b) where solid curves denote the result of the AMM 
[Eqs.(36) and (38)] and where marks express those of DS calculated by using H[x) — x. 
With increasing K, p is increased while 7 is little changed except for A^ < 5. 

3.1.2 Noise-strength (a) dependence 

The a dependences of stationary 7 and p for diffusive couphngs with A" = 10, A = 1.0 
and (3 — 1.0 are shown in Figs. 2(a), where filled and open marks denote 7 and p, 
respectively, in DS, and sohd and chain curves the respective results in the AMM. Note 
that the results of p are multiphed by a factor of 10 (= A^), and that three curves in the 
AMM are degenerated in Fig. 2(a): 10p(J = 0.5) = 10p(J = 0.0) = 7(J = 0.0). For 
J = 0.0, the relation of p = 7/A" holds in both the AMM and DS. For J = 0.5, 7 in the 
AMM is decreased compared to that for J = 0.0, in agreement with the result of DS. In 
contrast, p in the AMM is the same as that of J = 0.0 because it is independent of J [Eq. 
(25)], while p in DS is decreased with increasing J at o; > 0.6. The results of the AMM 
diverge at a = 1, where those of DS remain finite. This difference in p between the AMM 
and DS at large a is attributed to the second-moment approximation in the AMM, and 
it is the fallacy in the AMM neglecting higher-order moments. 

Figure 2(b) shows the a-dependent 7 and p for the sigmoid couplings of K — 0.5. 
Solid and chain curves denote 7 and p, respectively, in the AMM given by Eqs. (36) and 
(38): filled and open squares express those in DS with H{x) — x. Although 7 and p 
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diverge at a = 0.71 in the AMM, they persist up to a ~ 0.77 in DS. 
3.1.3 Coupling {J, K) dependence 

Figure 3(a) shows the J dependences of 7 and p for the diffusive couphngs with N = 10, 
a — 0.5 and (3 — 1.0: sohd and chain curves express 7 and p, respectively, in the AMM, 
and filled and open circle the respective results in DS. With increasing J, 7 is decreased 
while p is almost independent of J. 

The K dependences of 7 and p for the sigmoid couplings are plotted in Fig. 3(b) 
where solid and chain curves express 7 and p, respectively, in the AMM, and filled and 
open circle denote 7 and p, respectively, in DS. In the AMM, the critical coupling Kf. 
where 7 and p diverge, is = 0.75 [Eqs.(36) and (38)] while DS leads to ~ 0.84. 

These differences realized in numerical results are consistent with the items (ii) and 
(iii) mentioned in Section 2.3. We note in Figs. 3(a) and 3(b) that the synchronization is 
increased with increasing J and K because Eq. (41) approximately yields S oc (lOp — 7): 
S is proportional to the difference between chain (lOp) and solid curves (7). This agrees 
with the item (iv) in Section 2.3 denoting the similarity between the two couplings. 

3.2 Dynamical property 

We apply a pulse input given by 

&\t) = Ae{t-h)e{t2-t), (46) 

with A — 0.5, ti — 40 and ^2 — 50, where Q{x) denotes the Heaviside function: Q{x) = 1 
for a: > and zero otherwise. Figures 4(a), 4(b) and 4(c) show the responses of p(t), 7(t) 
and p{t), respectively, to the external input given by Eq. (46) with a = 0.5, /3 = 1.0, 
J = 0.0 and = 10. Solid curves denote the results of the AMM which are in good 
agreement with those of DS. Input pulse induces changes not only in p(t) but also in 7(t) 
and p{t). These arise from terms of p^ in Eqs. (24) and (25). Indeed, in the case of 
a = 0.0, input pulse induces no changes in ■j{t) and p{t), as shown by chain curves for 
the AMM result. 

The response of p{t) for the diffusive coupling is independent of the coupling J, as 
realized in Eq. (23). In contrast, the response p,{t) for the sigmoid coupling shows much 
variety depending on its coupling K. Figure 5(a) and 5(b) show p{t) and S{t), respectively, 
for various values of K with a — 0.5, (3 — 0.0 and = 10, when the pulse input given by 
Eq. (46) is applied. With increasing K, magnitudes of p{t) are increased at 40 < i < 50 
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during which the input pulse is apphed. It is interesting that the synchronization S{t) 
is decreased at 40 < t < 50 by an apphed input pulse which reduces the ratio of p/7, 
and then S{t) is much increased at t > 50. For K = 0.0, S{t) vanishes because p = 'y/N 
in Eq. (41). Figure 5(c) and 5(d) show similar plots of /i(t) and S{t), respectively, for 
combined noises of a = 0.5 and P = 1.0. With increasing K, the magnitude of fi(t) is 
again increased, although an agreement between the results of the AMM and DSs become 
worse than that shown in Fig. 5(a). S{t) is decreased by an applied pulse, but no increases 
at if: > 50, in contrast with the case shown in Fig. 5(b). 
We have applied also the sinusoidal input given by 



where A = 0.5 and Tp = 20. The responses of /i(t), 7(t) and p{t) are shown in Figs. 
6(a), 6(b) and 6(c), respectively, when the external input given by Eq. (46) is applied for 
a = 0.5, (3 = 1.0, J = 0.0 and N = 10. Solid curves expressing the results of the AMM 
are in good agreement with dashed curves of those of DS. Input pulse induces changes in 
fi{t) and also in 7(t) and p{t). For a comparison, we show, by chain curves, the AMM 
result for a = 0.0, for which no changes in 7(t) and p{t) by an applied input. 



It is interesting to discuss the stationary distribution of our generalized Langevin model 
given by Eqs. (l)-(3). In the case of no couplings {J = K = 0), the probability distribu- 
tion p{{xi},t) is given by 



where p{xi,t) expresses the distribution for a local variable Xi satisfying the FPE given 




(47) 



4 DISCUSSION AND CONCLUSION 



p{{Xi},t) = Ilip{Xi,t), 



(48) 



by 



-p(x„t) = - — {[F{x.) + ^G'ixi)G{x.) + &^]p{x.,t)} 

+ l^{[a'Gixif + p']pix.,t)}. (49) 



For a constant input of /*^^)(t) = /, the stationary distribution p{xi) is expressed by 

m-m 




(50) 
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with 



X{x) 
Yix) 



2 j dx 
2 f dr 



Fix) 



a^G{xy + /32 
/ 



(51) 
(52) 



_a^G{x)^ + (3^ 

For the hnear Langevin model with F{x) = —Xx and G{x) = x, p{x) in the Stra- 



tonivich representation becomes 

p{x) oc 



a 



-{X/aHl/2) 



exp[y(x)], 



with 



, (2I\ fax 
y[x)= — arctan — 
\afJI \ fJ ^ 



(53) 



(54) 



We examine the some hmiting cases of Eq. (53) as follows. 

(a) Equation (53) in the case of / = Y{x) = expresses the g-Gaussian 011111231211, 
which becomes, in the limit of large x (^ P/<y), 



p{x) oc x~ 



with 



S=— + l. 



The expectation value of x is given by 



/3' 



2{\-a^y 

which requires < A for positive definite (a;^). 
(b) For a = and 7^ 0, we get from Eq. (53) 



(55) 



(56) 



(57) 



p{x) oc exp 

(c) For P = and a 7^ 0, Eq. (53) becomes 

p{x) oc exp 



X — — 



A 







. {ah 


x_ 



{51 



(59) 



Figures 7(a)-7(c) show the distribution p{x) calculated with the use of Eqs. (53)- (59). 
The distribution p{x) for a = 0.0 in Fig. 7(a) shows the Gaussian distribution given by 
Eq. (58) without multiplicative noises, which is shifted by an applied input I. When 



12 



multiplicative noises are added {a ^ 0), the form of p{x) is changed but the average of 
(x) is conserved as shown in Fig. 7(a). Figure 7(b) shows that when the magnitude of 
additive noises (3 is increased, the width of p{x) is increased. We note in Fig. 7(c) that 
although p{x) is symmetric for 1 = 0, the external input / increases the asymmetry in 
p{x). Figures 7(a)-7(c) clearly show that p{x) is much modified by the presence of /. 

Now we consider the averaged, global variable X{t) given by Eq. (8). The stationary 
distribution for a global variable X given by Eq. (9), is analytically expressed only for 
limited cases. 

(a) For /9 7^ and a = 0, P{X) is given by 

fxrr 



P{X) oc exp 



X 



A 



which arises from the central-limit theorem for 
(b) For / = 0, we get 



P{X) = — dk e'^^ <^{k), 
2tt J-oo 



with 



$(fc) 



N 



n N 



where 0(/c) is the characteristic function for p{x) given by 



m 



-ikx 



p{x)dx, 



K^{\' I k 



with 



(60) 



(61) 
(62) 

(63) 
(64) 



A' 



a 



(65) 
(66) 



Ky{x) expressing the modified Bessel function. 

The asymptotic form of P{X) for large X and large N is obtained as follows. By using 
the relation: 



z'^B^{z) oc 



cz 



dz^" + 



with 



1 



for 2; -C 1, 1/ 7^ 1 



d 



%-!)' 

^.2v Yd 



V) 



2 J r{i + u)' 



(67) 

(68) 
(69) 
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we get, for large A^, 

$(/c) cx exp(-aiv k"^), for z/ > 1 (70) 



with 



oc exp(-6^ I k r), for < z/ < 1 (71) 



= cN-\\'f, (72) 
bj^ = dN^-^^iX')^". (73) 



For large X, Eqs. (61), (70)-(73) yield 



X2 

oc X-^\ for < z/ < 1 (75) 



P{X) oc exp --^ I , for z/ > 1 (74) 

\ ^0" A 



with 



5' = 2z. + l = 4 + l- (77) 

It is interesting that for = 1, Eq. (76) coincides with Eq. (57) and the index of 6' 
given by Eq. (77) is the same as 6 given by Eq. (56). The case of z/ = 1, excluded in the 
above analysis, will be numerically studied below. The stable distribution of P{X') for 
X'{t) = Y.i Xi{t) with < z/ < 1 was discussed in Ref. |^. 

Figures 8(a) shows distributions of a global variable P{X) for I = 0.0 calculated by 
DS for the Langevin model given by Eq. (1) with = 1, = 10 and A^ = 100 (A = 1.0, 
a = 1.0, P = 0.5 and u = 1.0). For a comparison, results of the analytic expression 
given by Eqs. (61), (62) and (64) are plotted with a shift by X = —2 for a clarity of the 
figure. We note that with increasing A^, the width of P{X) becomes narrower, which is 
consistent with the central-limit theorem. Figures 8(b) and 8(c) show P{X) for / = 1.0 
and I = 2.0, respectively, calculated by DS. The A^ dependence of P{X) for finite I is 
intrigue: with increasing A^, not only its width becomes narrower but also its profile is 
considerably modified, as shown in Figs. 8(b) and 8(c). This trend is more significant for 
a larger /. 

So far we have assumed the vanishing couplings, which is now introduced. Figure 
9(a) shows distributions of p{x) and P{X) for the diffusive couplings of J = 0.0 (dashed 
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curves) and J = 1.0 (solid curves) with / = 0.0. We note that with increasing J, the 
width of p{x) becomes narrower while that of P{X) is not changed. This is the case also 
for finite / = 1.0, as shown in Fig. 9(b). 

Figure 9(c) shows p{x) and P{X) with / = for the sigmoid coupling of K = 0.0 
(dashed curve) and K = 0.5 (solid curve). We note that the width of P{X) for K = 0.5 
become wider than that for K = 0.0. Figure 9(d) shows that an introduction of K with 
finite / = 1.0 induces not only an increase in the width of P{X) but also its shift. This 
is in contrast with the case of the diffusive coupling shown in Fig. 9(b), where P{X) has 
little effects of J. 

The coupling dependences of stationary distributions of p{x) and P{X) are related to 
those of 7 and p, because ^7 and y^p approximately express the widths of p{x) and P{X), 
respectively. Figure 9(a) shows that with increasing J, the width of p{x) is decreased 
while that of P{X) is unchanged for the diffusive couplings. In contrast, for the sigmoid 
coupling, the width of P{X) is increased while that of p{x) is unchanged when K is 
increased, as shown in Fig. 9(c). These are consistent with the dependences of 7 and 
p on the type of couplings expressed in the items (ii) and (iii) having been discussed in 
Section 2.3. 

5 CONCLUSION 

By using the AMM, we have studied stationary and dynamical properties of finite A^- 
unit Langevin model which is subjected to multiplicative noises and which is coupled 
by diffusive and sigmoid couplings. Properties of coupled Langevin model are shown 
to depend on both the type and magnitude of couplings, which is the same as in the 
case of FitzHugh-Nagumo model jU] [21] • One of advantages of the AMM is that we may 
easily solve the low-dimensional differential equations although its applicability is limited 
to the weak-noise case. It would be interesting to apply the AMM to various types of 
stochastic coupled ensembles such as neuronal and complex networks in order to discuss 
their dynamics. 
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Figure 1: (Color online) (a) The N dependences of the stationary 7 and p for various 
diffusive couplings (DC): circles, triangles and square denote results of DS for J = 0.0, 
J = 0.2 and J = 0.5, respectively, and solid curves express those of the AMM. (b) The N 
dependence of the stationary 7 and p for various sigmoid couplings (SC) with H{x) = x: 
circles, triangles and square denote results of DS for K — 0.0, K — 0.2 and K — 0.5, 
respectively, sohd curves express those of the AMM (A = 1.0, a — 0.5 and /3 = 1.0). 
Dashed curves are drawn only for a guide of the eye. 

Figure 2: (Color online) (a) The a dependences of the stationary 7 (solid curves) and p 
(chain curves) for the diffusive couplings (DC): circles, triangles and square denote results 
of DS for J = 0.0 and J = 0.5, respectively, and solid and chain curves express those of 
the AMM. (b) The a dependence of the stationary 7 (solid curves) and p (chain curves) 
for sigmoid couplings (SC) oi K — 0.5 with H{x) — x: squares denote results of DS, and 
solid and chain curves express those of the AMM. (N = 10, A = 1.0 and P — 1.0). Dashed 
curves are drawn only for a guide of the eye. 

Figure 3: (Color online) (a) The J dependences of the stationary 7 (solid curves) and p 
(chain curves) for the diffusive couplings (DC): circles denote results of DS, and solid and 
chain curves express those of the AMM. (b) The K dependence of the stationary 7 (solid 
curves) and p (chain curves) for sigmoid couplings (SC) with H{x) — x: squares denote 
results of DS, and solid and chain curves express those of the AMM. {N = 10, A = 1.0, 
a — 0.5 and P — 1.0). Dashed curves are drawn only for a guide of the eye. 

Figure 4: (Color online) Responses of (a) p{t), (b) 7(t) and (c) p{t) to the pulse input 
with J = 0.0, = 10, A = 1.0, and jS = 1.0: solid and dashed curve denote results 
of AMM and DS, respectively, for a = 0.5, and chain curves express that of AMM for 
a = 0.0. 

Figure 5: (Color onhne) (a) Responses of p{t) and (b) S{t) to the pulse input for various 
sigmoid couphngs with A = 1.0, a — 0.5, /3 = 0.0 and N — 10. (c) Responses of p{t) and 
(d) S(t) to the pulse input for various sigmoid couplings with A = 1.0, a = 0.5, /? ~ 1.0 
and N — 10. Solid and dashed curves denote results of AMM and DS, respectively. 

Figure 6: (Color online) Responses of (a) p{t), (b) 7(t) and (c) p(t) to the sinusoidal 
input with J = 0.0, = 10, A = 1.0, and [3 = 1.0: solid and dashed curve denote results 
of AMM and DS, respectively, for a = 0.5, and chain curves express that of AMM for 
a = 0.0. 

Figure 7: (a) Distributions p{x) of local variable x for various a with A = 1.0, P — 1.0 
and / = 1.0, (b) p{x) for various (3 with A = 1.0, a — 1.0 and / = 1.0, and (c) p{x) for 
various / with A — 1.0, a — 1.0 and P — 0.5. 
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Figure 8: (Color online) Distributions P{X) of global variable X calculated by direct 
simulation (DS) for = 1 (dashed curves), = 10 (solid curves) and = 100 (chain 
curves) with (a) 7 = 0, (b) / = 1.0 and (c) I = 2.0 (A = 1.0, a = 1.0 and P = 0.5). 
Results calculated with the use of Eqs. (67), (68) and (70) are plotted in (a) with a shift 
oi X — —2 for a clarity of the figure. 



Figure 9: (Color online) Distributions of p{x) and P{X) of local and global variables, 
respectively, with (a) / = 0.0 and (b) / = 1.0 for diffusive coupling (DC) with J = 0.0 
(dashed curves) and J = 1.0 (solid curves), and those with (c) / = 0.0 and / = 1.0 
for sigmoid coupling (SC) with K = 0.0 (dashed curves) and K = 0.5 (solid curves). 
{N = 10, A = 1.0, a = 0.5 and P = 1.0). 
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